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ABSTRACT 


The idea of survivability introduced as a network parameter has led 
to so many investigations. Several measures of survivability has been 
studied. The number of links and/or stations needed to be damaged to 
disrupt the system is the survivability criterion adapted in this study. 

The development of analysis procedures for directed, undirected, or 
mixed networks based on the above criterion &ànd use of the concepts in 
network flow and graph theory are treated in detail including computer 
program implementation of the algorithms, Finally a practical design 
algorithm for minimum-cost survivable network with respect to branch dis- 


connection using a heuristic approach and analysis techniques isdescribed. 
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I. INTRODUCTION 


Network complexes such as railroads, highways, communication circuits, 
telephone systems, and many others which are interconnections of links and 
stations, and where there exist commodity flow, have been modelled by lin- 
ear graphs. The vertices of the graph symbolize stations or terminals 
while the branches denote lines, trunks, or circuits, The graphical struc- 
ture of such systems allows the use of topological concepts in the develop- 
ment of the theory of network. 

In 1961 the idea of survivability as a graph parameter was introduced. 
This concept was later adapted as a network parameter. Survivability of 
a network is defined as its capability to function as a system after dam- 
age induced by enemy attack, Survivability and vulnerability are used 
interchangeably by many authors. The term network reliability is more ap- 
propriately used when referring to damage caused by natural disturbances. 

The study of network survivability employs topological methods and the 
concept of network flow. There are two distinct areas, namely: determin- 
istic survivability and probabilistic survivability. A network is con- 
sidered probabilistic when the parameters involved like commodity flow and 
weights assigned. to links and stations are statistical. When there is 
certainty in the existence of these parameters, the network becomes de- 
terministic, 

Up to present, network survivability has been the subject of research 
by many graph and network theoreticians. Recent works were focused on the 
Mathematical formulation of survivability criteria, the development of 
analysis methods, and the synthesis of optimally invulnerable networks, 


Numerous results have already been published. 





This thesis is based on the paper by Frank and Frisch [10 | which 
summarizes the most significant results of the researches in this area. 
The treatment of deterministic networks is restated in more detail and ex- 
tended by adapting other results in network flow theory with some modifi- 
cations to provide additional analysis procedures. To provide working 
tools for analysis and synthesis, the existing algorithms and those that 
were modified were coded for the computers, Further simplifications were 
introduced when encoding these algorithms into the programs which appear 
in the appendices attached herewith. 

In Chapter II of this report, a summary of the important elementary 
principles of network theory essential in the study of network surviva- 
bility is presented, Chapter III covers the analysis procedures based on 
some survivability criteria, The principles from which they were derived 
are also discussed. Investigation of directed, undirected, and mixed net- 
works with respect to branch disconnection and vertex disconnection are 
separately treated. 

These analysis tools may be employed in the design of optimally sur- 
Vivable networks, The object of Chapter IV is design optimization by an 


heuristic approach using analysis techniques, 





II. NETWORK FUNDAMENTALS 


2.1 NETWORK TOPOLOGY 

One of the many applications of linear graphs to engineering problems 
is the topological representation of networks. Networks may be categorized 
as "electrical networks", "switching networks", and "flow networks", The 
latter which is often called "communication net" has associated with its 
branches some information or commodity flow. Symbols, terminology, E a 
minimum of the fundamental concepts necessary in the study of this type of 
network is discussed in this chapter. 
2.1.1 Directed, Undirected, and Mixed Networks 

A directed network is represented by a directed linear graph G-(V,I'), 
consisting of a set of elements called vertices denoted by the symbol V 
and a set of BE pairs of vertices called branches denoted by |. The 
vertex Ve is called the initial vertex and the vertex E the terminal ver- 
tex; vertices Vs and u are said to be adjacent. 

A directed branch is indicated by the ordered pair (i,j). Each 
branch carries an orientation and the information flow is only in the spec- 
ified direction, 

An undirected network is m by an undirected graph G=(V,ff), 
consisting of a set of vertices V, and a set of unordered pairs of these 
vertices denoted by ff. An undirected branch is indicated by ЕД, өзе 
E! The flow of information is bidirectional. 

A mixed network is symbolized by a graph that has both directed and 
undirected branches. A non-oriented branch may be represented graphically 


by two parallel branches oppositely directed. 





The network graph may have more than one branch between the same 
pair of vertices. 
2.1.2 Paths In a Network Graph 

À directed s-t path T. is a sequence of distinct vertices and 


branches connecting two arbitrary vertices, = and v It has the form 


t' 
e (5,1), p. (1,,1,), En e.e. у d (i, t), De 


Àn undirected s-i path T. is a sequence of distinct vertices and 
branches of the form 

vo L5). ], v, С Y v. 

| 2 k 

Two directed or undirected paths ТІ, ала МЕ are branch disjoint if 
they have no branches in common and vertex disjoint if the only vertices 
they have in common are Ја and Vie 

In a sequence of vertices and branches describing a directed path, 
if (i,j) is in the path, then (i,j) is a forward branch with respect to 
the path. A branch (j,i) in the path is a backward branch with respect 
to the path, 

The network graph G is said to be connected if for any pair of ver- 
tices Va and С there is a path Mos in G, Otherwise it is disconnected. 
The resulting Марне are dioan 
21,3 Cutsets and Cuts 

An ect ed branch cutset of an undirected, connected graph is 
the minimal set of branches, the removal of which yields a graph of two 
or more components, 

A directed branch cutset is a minimal set of branches of a directed 
graph, the removal of which disconnects all dürected paths from at least 


one vertex of G to at least another vertex of G. 





Àn undirected vertex cutset of an undirected, connected graph is 
the minimal set of vertices whose removal from G disconnects the graph. 

А directed vertex cutset is a set of vertices of a directed graph 
whose removal from G destroys all directed paths from at least one of the 
remaining vertex of G to at least another remaining vertex. 


Cutsets with respect to a vertex pair „з and v, are referred to as 


t 


Бор Cutsets, 


A cut 259 separating М2 ала Уһ is a set of branches (i,j) such 
that а € X and Ст Gar, In Riese cif ~ У ОМОТ and X = М 
VıV, | » the set of branches (X,X) = [n (1225 (2,3) is a directed 
s-t cut. Branch (2,4) is oppositely directed with respect to the vertices 


s and t therefore it does not belong to the set. 





Ш 
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2.1.4 Network Tree 

А spanning tree T = (V,U) of a connected network graph G = (У, ) 
is a connected subgraph of G containing all of the vertices of G but no 
closed path. 

If n is the number of vertices in G, a tree contain n-1 branches. 

Every pair of vertices vi and МЕ of graph T is joined by one and 
only one i-j path. Therefore it becomes disconnected when one branch is 


removed. 


2.2 FLOWS IN DETERMINISTIC NETWORK 
2.2.1 Branch Flow and Branch Capacity 

Aside from the topological structure of a communication net, there 
is assigned to each branch a finite weight f(i,j) or RE | called the 
branch flow from v, to Vee Also a weight c(i,j) or cli,j] is assigned to 


each branch which is the maximum amount of commodity that can reach а 


from v. per unit time, called the branch capacity. 
For an undirected network, 
c[1,3] » -cL3,1], (CENT) 
£[i,3] CBE NE (252210) 
a a e (2,2305) 


are any chosen source and sink vertices respectively, For 


|| 


where v and v 
25 t 
a directed case, equations (2.2.1.1) and (2.2.1.2) are.not necessarily 
true, 
An undirected branch can be represented by two oppositely directed 


branches between А апа E With capacities 


Cm B c(j,i) = CS iE (2 ен) 
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2.2 Flow Pattern 
The set of flows associated with the branches in G is called a flow 
pattern denoted by F, A flow pattern is said to be feasible if it satis- 


fies Equations (2.2.2.1) and (2.2.2.2) for some nonnegative constant T 


ie 
3 
which is called the value of F. 
For all Vs € V of a directed network, 
mE а= 
ANS у =н МЕШИ Е 2282.1) 
FO i= 
О < £(4,5) < eG ay тота eee (2027202) 


This means that the net flow out of the source А 18 T A and the 
) 


flow out of the sink vs is E ; whereas the net flow out of an intermed- 


yt 


late vertex is zero, This latter relation is also true for an undirected 


network, 


|| 


A branch (i,j) is saturated if f(i,j) » c(i,j) and unsaturated 


0. 


|| 


Otherwise. It is said to be empty if fli,j) 
2.2.3 Residual Capacity and Augmentation Path 

The residual capacity (Tg 4) is defined both for directed and un- 
directed network as 


em, 4) - Кіп Е [ £(3,3) » Mim 4 c(a,b) - ау] | 


backward forward 
branches branches 
Қылау S "Sit (a,b) € Wr C202. 3 8) 


A path 1, , is saturated if r(m. 2 = 0 and unsaturated otherwise, 
) 


t 
, 
Additional flow can be sent through an unsaturated path obtaining a new 
mOn pattern. To achieve this, add r(m. D to all forward flows and 
) 
subtract r(m. D) from all backward flows along the path. 
3 


If à path between two vertices has 


Uc 
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and fho 
then we call this a flow augmentation path. The flow along this path can 
be increased according to the previously stated rule. 
2.2.4 Maximizing Network Flow 
One of the most important and basic theorems in network flows is 
the "Max-flow Min-cut theorem" formulated by Ford and Fulkerson Ди The 
problem of obtaining the maximum obtainable flow value - through a net 
with reference to two vertices у. and v, is solved using this theorem. 
Max-flow Min-cut Theorem: The maximal flow value obtainable in a 
network G form Ya to E is the minimum of the value of a cutset taken over 
all cutsets separating va and Vie 
To state the theorem mathematically, let: 


" t be the value of the minimum s-t cut 

j 

с(Х,,Х,) be the capacity of a cut (X,,X.) 
EI vr 

Год) | be the set of all s-t cuts in the net 
0-17 вур 


then: 


= Een un Le(X, X, ), +] = е” (ZWAR, 


The value v of the cut (USES) is defined as 
у = f(X,X) - £(X,X) « e(X,X) (Bea) 
The Ford and Fulkerson Labeling and Augmentation algorithm [7| based 


on the above theorem provides a method of increasing f systematically 


ar 


from an arbitrary feasible flow pattern with value less than gra ү and ter- 
5 
Minates when f С а T. 4 The resulting flow pattern is maximal. 
U j 


The algorithm will be stated in a modified form in Chapter III. 
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2.3 MATRICES ASSOCIATED WITH SURVIVABLE NETWORKS 
| Network graphs may be structurally drawn or be described by matrices. 
The latter representation is useful for computer calculations. This sec- 
tion enumerates several matrices used in the study of survivable networks. 
2.9.1 The Incidence Matrix 
Let G — (V, f! ) be the directed graph of the network consisting of 
n vertices and m branches. Let A = (а, 4 be ann x m matrix whose i-jth 
entry is а, D Then A is the incidence matrix of G if 
+1 if a branch is directed away from V; 
а. . = 4-1 if a branch is directed towards Vi 
O otherwise 
If G is an undirected graph, A = (а, 4 is the incidence matrix if 


1 GP oratith is in@ident av Vs 


0 otherwise 
2 Connection Matrix 
An alternate representation of a non-oriented network is the con- 
nection matrix K = Lk; 3) K is ann xn symmetric matrix whose i-jth 
3 
entry is k, . where 
EL 1j 


ИЕ Сап 


O otherwise 
For a directed network, a -1 is entered for an oppositely directed branch. 
2.3.3 Branch-capacity and Terminal Capacity Matrices 
The other way of representing a communication net of n vertices is 
by using a symmetric matrix C = іс, 4 called the branch-capacity matrix 
Where с. . = c. , is the sum of the capacities of the branches connected 


у Jat 


between Vi and л For an oriented network, the C-matrix is asymmetrical. 
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For any two vertices in G, there exist a finite, maximum obtainable 
capacity of transferring information from one vertex to the other called 
the terminal capacity. The terminal capacity matrix T = [ts 4 represents 
the maximum flow that can exist between all the possible vertex pairs of 
G. The T-matrix is also asymmetrical for a directed network, 

2.3.4 Connectivity and Redundancy Matrices 

The connectivity matrix is essentially a T-matrix with each branch 
of G or G assigned a maximum capacity of unity. 

Ann xn matrix R = EM is a redundancy matrix where the i-jth 
entry ias is the number of branch/vertex disjoint paths between V; and 
И It is symmetrical for a non-oriented net, А uniform R-matrix is one 
which has identical entries except the main diagonal. 

The main diagonals of these matrices have infinity as entries how- 


ever for computer calculation zero entries are used, 


2.4 DISCONNECTING SETS OF A NETWORK 

A connected network represented by G can be disconnected by deleting 
certain branches or vertices, 

An s-t branch cutset is a set of branches whose deletion breaks all 


directed paths from um to v,; likewise an s-t vertex cutset is a set of 


49 


vertices whose deletion breaks all directed paths from ve to v A mixed 


+ 
s-t cutset is composed of branches and vertices, 

Clearly if the capacity of each branch or vertex is set to unity, 
then the value of the cutset is numerically equal to the number of com- 
ponents of the cutset, 

Let Wt be the number of elements in the smallest s-t vertex cut- 


set; 


Т. t be the number of elements in the smallest s-t branch cutset; 
, 
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т. ¿ be the number of elements in the smallest s-t mixed cutset; 
j 


then the disconnecting sets of the connected graph are as follows: 


W) = Min LW. +] (ДЕШІ) 


S,UC V 
E 


T = Min Кы (Omm 


St C V 
sft 


Т = Min Kes xi | (ОСИК ) 


st E Vy 
sft 


The smallest disconnecting set of à connected graph denoted by v is 


defined as 


IA 1 (2.4.1.4) 
However it was shown in [3] that 
у == Мап [| 6) m le) 
cL CV 
SD 
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IIl. ANALYSIS OF DETERMINISTIC SURVIVABLE NETWORKS 


В.1 SURVIVABILITY CRITERIA 

A significant criterion in network analysis is its survivability or 
vulnerability. Іп а communication net, this parameter refers to the de- 
gree in which the system remain functional after an enemy attack or nat- 
ural disaster. 

No general survivable criteria has yet been formulated however several 
measures of survivability has been investigated namely: 

1. The number of links and/or stations that must be destroyed before 
communication is disrupted [4], 

2. the minimum number of links that must be deleted in order to iso- 
late a group of stations from communication with the remainder of the net 
[4]; 

3. the minimum number of isolated subnets that would result from an 
attack or destruction. This is termed the "independence number" of the 
network [1], 

4, fraction of stations which survive an attack [2], [6], 

5. fraction of stations that can be reached from a given point by 
a path of no more than r links after an attack [8]. 

The first three criteria apply to deterministic network while the 
last two criteria are associated with probabilistic net. 

Analysis of large networks based on the first criterion is treated 
in this thesis. This index of survivability suggests the application of 
the concept of disconnecting sets. More specifically, the parameter (y 
is a measure of vulnerability to station destruction and 7 is a measure 


to link destruction. 


Le | ы 





3.2 ANALYSIS WITH RESPECT TO BRANCH DISCONNECTION 
3.2.1 s-t Branch Connectivity 

The problen discussed here is the finding of the maximum number 
of branch disjoint s-t paths between any given pair of vertices ya and 
1% in G. Setting the branch capacities to unity and applying the max- 
is obtained 


flow min-cut theorem, the maximal flow f from n to v 


n t 
which is numerically equal to the smallest branch cutset А This 
value further represents the maximum branch E paths between the 
given vertex-pair. 

The Ford-Fulkerson algorithm solves a general max-flow problem be- 
tween two vertices, the source and sink. It is stated here in simplified 
form to adapt it to the special case being investigated. 

A few terms are defined below before stating the algorithm: 

1. A vertex is labeled by an ordered pair (i,+). For a particu- 
lar vertex ОЛ the index 1 denotes the preceeding labeled vertex E 
The second term in the label is plus if for a forward branch, f(i,j) = 0 
and a minus sign if for a backward branch, f(j,i) > 0. Otherwise a ver- 
tex is unlabeled. 

2. A paretroulbssavernbtes ш, is scanned if its adjacent vertices E 
are labeled from it, Otherwise ya is unscamned, 

З. А vertex is considered to be in one of the three states, namely: 

a. unlabeled and unscanned 
b. labeled and unscanned 
с. labeled and scanned 

The following steps determines the s-t branch connectivity of a 
directed or mixed network: 

1. Represent the network by a directed linear graph. Number the 


vertices consecutively from 1 to n. 
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2. For each branch (i,j) € FF, assign c(1,3) = 1 and 21,3) een 
3. Choose the source and sink vertices ue and Ye 
4, Apply Algorithm 1 
Algorithm 1 (Modified Ford and Fulkerson Algorithm) 
Labeling Routine: 
a. Erase all vertex labels 
b. Label E by (0,+). Source is now labeled and unscanned, 
all other vertices are unlabeled and ace 
с. Select any labeled unscanned vertex v, (including у) 
Assign labels to all unlabeled vertices Е adjacent to У. such that 

i) for forward branch (i,j) with flow f(i,j) = 0, label 
Va by (1,7). Vertex v; is now labeled and unscanned, 

2) for backward branch (j,i) with flow f(j,i) » O, label 
v; lye (i,-). Vertex У. is now labeled and unscanned. 

With all adjacent vertices to Vi labeled and unscanned, 

V; is now scanned, 
а. Repeat Step c until either 

1) Vy is labeled or 

2) no more labels can be assigned and V4 is unlabeled, 

In Step a(1), an augmentation path is found and an increase 
of flow is possible. Proceed to Augmentation Routine. 

In Step d(2), the algorithm terminates. Flow is now max- 
imized. The disconnecting cutset are the branches joining the labeled 
and unlabeled vertices. The value of the cut is the branch connectivity. 

Augmentation Routine: 
a. Let z = t and go to Step b. 


b. If label on v, is (q,*), set f(q,z) = 1. If label on УЕ 


Eu set #(2,а) = 0, 
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Ce If q = s return to Step a of the Labeling Routine. Otherwise 
set z = q and go to Step b of the Augmentation Routine. 

Example: A mixed network 

Fig. 3.1(a) is the linear graph representation of the network, 
In (b), the non-oriented branch (3,6) is split into two oppositely direct- 
ed branches (3,6) and (6,3). The first labeling obtained the augmentation 
Patch s,(s,3),3,(3,7) ,7,(7,t). In (c) the second labeling found the aug- 
mentation path s,(s,),4,(4,3),3,(3,6) ,6,(6,t). In (d) the sink vertex 
v, cannot be labeled since (6,8) ana (7,8) are already saturated and (5,8) 
is a backward branch with zero flow. The algorithm terminates. The cut 
is {(6,8),(7,8) | and s-t connectivity is two. 

Flow chart for this analysis is in Appendix A-1 and the computer 
program MAXCON-1A is in Appendix B-2. То minimize computer storage, a 
vertex label is further modified by using a single signed integer (+1) 
instead of the ordered pair (i,+). When i is negative the magnitude is 
obtained by just multiplying by minus i. 

In the program, vertex scanning is accomplished by using two stor- 
age vectors. The first vector is a list of the labeled vertices being 
scanned while the second provides storage for the vertices being labeled. 
Once the scanning process of all the vertices listed in the first vector 
is completed, the contents of the second is transferred to the first which 
in turn becomes the list of vertices to be scanned, 

Algorithm 1 is applied to the analysis of undirected network with 
slight modification. For a non-oriented net, each branch flow is con- 
sidered bidirectional with a maximum capacity of unity. However, in no 
case can a positive back-flow occur since each augmentation path has as 


its elements all forward directed branches with respect to A and Vie 
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Figa а! 
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Consequently there cannot be any negatively labeled vertex. Therefore the 
tests for a backward branch and а negative-labeled vertex, x adjust- 
ment of backward flow are not necessary in this case. 

These steps are deleted in Algorithm 1 and MAXCON-2A (Appendix 
B-4) is developed for determining s-t branch connectivity of an undirected 
network. 

3.2.2 Multiterminal Branch Connectivity 

For an n-station network, its graph "m have rS possible 
vertex pairs. To apply Algorithm 1 n(n-1)/2 times in order to obtain the 
connectivity between all the possible vertex-pairs would be lengthy. 

Gomory and Hu [13] have shown that for an undirected network with 
n vertices, any flow in the network is numerically equal to some flow in 
a maximal spanning tree, Since there are only n-1 branches in a spanning 
tree, then there can only be n-1 numerically different flows possible. 
Therefore all the n(n-1)/2 possible maximal flows that can exist in the 
network can be deduced after doing n-1 flow maximizations. 

Applying the above principle and the max-flow min-cut theorem in 
the analysis of non-oriented nets, the multi-terminal network connectivity 
problem can be solved by only n-1 separate applications of Algorithm 1. 
The analysis procedure is embodied in the following algorithm. 

Algorithm 2 

1. Choose two vertices Уз апа Vy in the undirected graph G. 
With branch capacities set to unity, find an initial s-t branch cut using 
Algorithm 1. Represent this cut by a generalized vertex tree with con- 
densed vertices X and X connected by a branch with capacity equal to the 
computed cut value бе: 


2, Draw a condensed network graph G, where X is represented 


1 


by a single vertex. Select two vertices V and Ya in X and find the 
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smallest К-1 cut by Algorithm 1. Form the generalized vertex tree by 
dividing X into X, and X, with 7°, as the capacity of the branch con- 


necting them. In this tree, X is adjacent to X, if it is in the same 


1 


side of the cut as X,, or X is adjacent to X, if it is in the same side 


of the cut as X,. 
3. The above procedure is repeated. With each step, a tree 
of generalized vertices connected by branches with corresponding cut 
values is obtained. To proceed with the computation, 
ae Select a generalized vertex X. 
b. Condense each component in G except X. into a single 
vertex forming the condensed network graph Ge 


с. Choose two original vertices 25 ала у, іп X. as source 


i 
and sink and apply Algorithm 1 to с, to get Er T 
The cut obtained in this process splits X. into two parts 

Е, апа E This is represented in the generalized vertex tree by re- 
placing X. by the condensed vertices = and E. connected by a branch 
bearing the cut value, All other branches and vertices in the tree are 
unchanged except those which were formerly connected to X. Apply the 
adjacency rule as in Step 2, 

l, Repeat the process until these generalized vertices con- 
sist of E one vertex each, This results in a spanning tree. Thus 


for any chosen ШЕ and у,, the s-t connectivity is the smallest branch 


12 
weight in the unique s-t path in the tree, 

An example for the application of the above procedure is given 
in Fig. 3.2. The generalized vertex tree and the condensed network с, 


resulting in each flow calculation are indicated. For the 12-vertex net- 


work, the spanning tree is formed after 11 max-flow computations. 
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n ivewchast for this analysis procedure is given in Appendix А-2, 
This is further encoded into a computer program entitled MAXCON-2 appear- 
ing in Appendix B-3. In this program the printout is not a spanning tree 
but rather a symmetrical branch-connectivity matrix which is essentially 
a terminal capacity matrix. Its entries are derived from the cut values 
appearing as branch weights of the trees generated after every max-flow 
calculation where a cut isolates a single vertex from the rest of the con- 
densed network graph G,. The connectivity BEE this single generalized 
vertex and all the rest of the vertices in G are the updating entries. 

Other information relevant to MAXCON-2 are stated in the program 
documentation. 

The concept of spanning tree cannot be applied to the analysis of 
a directed or mixed network as generally these networks are not symmet- 
trical, To obtain a multiterminal branch-connectivity matrix for an ori— 
ented net, it is necessary to do S flow maximizations by Algorithm 
1. Up to present, there is no avallable method that could reduce this 
number of calculations however Frisch has developed the Flow Variation 
Algorithm [12] which effectively reduces the number of steps required to 
maximize E 

- The ideas underlying this new algorithm are as follows: let Vos 

Vo and V4 be three arbitrary vertices in a graph G. Rather than apply 
Algorithm 1 to independently maximize f and f 


5% st у 


by Algorithm 1 then maximize Lo t in the usual manner until a vertex v. 
) 


is labeled for which f(g,V) > 0. Then use the flow pattern D(s',t) to 


first maximize f 


па F(s,t). 
Applying this basic idea to solve the multiterminal connectivity 
problem which is a special case of a general flow problem, all that is 


necessary is to slightly modify the algorithm to fit into the problem. 


P 





The result is a simplified version of the original algorithm. The following 


notations are used: 


D - flow path for vertex-pair Um and in 

et - flow path for vertex-pair Vor and 42 
F(s,t) - flow pattern for vertex-pair v, and v, 
D(s',t)- flow pattern for vertex-pair Mn and Ше 


f(i,j) - a branch flow in path m. 


nt 
d(i,j) - a branch flow in path pm 
9 
A 
ТІ - section of path т from v to v called "Augmenta- 
5,6 S,t S g 
irons pn 
пі, = section of path 7 , , from v-, їо у саше ите 
S',g s',t S g 
tion Раљћ" 
H : - . 
D - section of path Е from Уш to un called "Homing 
Path" 
І" m а loop detected by the Homing Routine; m is the first 
) 
vertex in the loop 
і, р 7 a loop detected by the Truncation Routine; p is the 
, . 


first vertex in the loop 
The subscript c whenever it appears with the symbols indicates the 
2 stage of max-flow calculation. 
The following routines are employed in the algorithm: 
Routine A - Detecting an Augmentation Path 
1. Label vI by (s,+). ve is now labeled and unscanned, А11 
other vertices are unlabeled and unscanned. 
2. Select any labeled and unscanned vertex Vis Assign labels 


to all unlabeled vertices adjacent to Vi such that 
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label ү. by Û, 
label У; С 


is now scanned, 


3. 


a. for forward branch (i,j) € T with flow f(i,j) = 0, 
+), Е is now labeled and unscanned, 
b. for backward"branch (j,i)-€ [98 mith flow (j i) > ӨҢ 
E v; is now labeled and unscanned, 


With all adjacent vertices to Va labeled and unscanned, Vi 


Repeat Step 2 until either: 
а. V. is labeled or 
b. a vertex - is labeled with à (g,V) > 0 


с, no more vertices can be labeled and vg is unlabeled, 


Observe that except for Step 3(b) this routine is the Ford 


and Fulkerson Labeling Routine, 


Routine 


1. 


2. 


1 - Detecting a Truncation Path 
Letz-g 


Find one vertex Va such that 4 (1,2) > 0 and label Va by 


еру ===> 
a, ІТ МЕ has assumed a value i before, terminate. 
b, ЈЕ = s', terminate. 


c. Otherwise return to Step 2. 


If the routine terminates in Step 3(a), a loop г. : has been 
3 


detected rather than a path. 


Routine 


1. 
2, 


(а) 


H - Detecting a Homing Path 
Letq=g 


Find one vertex ү; such that d(q,j) > 0 and label у; by 


Let gq = J 


a. If q has assumed the value j before, terminate. 
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b, If j = t, terminate. 
с. Otherwise return to Step 2. 
If the routine terminates in Step 3(a), a loop Ln has been 
detected instead of a path. 

The main algorithm for obtaining the multiterminal branch-connec- 
tivity of an oriented network is now stated in terms of Routines A, T, 
and H. 

Algorithm 3 

1. Choose two vertices a and ME Using Routine A, find an 
initial s-t cut, 
2. Using Routine A find an s-g augmentation path - or an 


S-t augmentation path m In the former cas® go to Step 3, Min the dare 


+" 
3 
ter case go to Step 8. If no such path can be found, terminate, 


3. Using Routine H find a Homing Path m 


a or a Homing Loop 
3 


г bs In the former Case go to Step 4; in the latter case go to Step 
j 


5(a). 


l|, Using Routine T find a Truncation Path m. or a Trunca- 


EU 
tion Loop boy: In the former Case go to Step 6; in the latter case go 
to Step 5(b). 

5, a. Set to zero all à (1,3) for all branches (i,j) in the 
loop m Erase all vertex labels and go to 2. 

b. Set to zero all à (1,3) for all branches (i,j) in the 

loop E. Erase all vertex labels and go to 2. 

6, Trace and erase g-t path, For every branch (en те 


set 
£ 40:3) = 1 


er) 0 


20 | : 





7. Truncate the s'-g path. For every branch (i,j) in п: 


, 


",Е 
сез НЕЗ = 0, 


8. Augment the s-g path, For every branch (i,j) in LM set 
£453) ce 
Erase vertex labels and return to Step 2. 
If a "loop of flow" is detected, the flow in this loop must be 
set to zero as in Step 5 as otherwise a nonfeasible flow will be achieved. 
As an example, consider the graph in Fig. 3.3. Let the indicated 


flow pattern D(s',t) be the result of initial maximization of 1 
9 


Designate у ав the source for the second stage of maximization. Step 2 


of the algorithm labels Ке and detects it as a vertex E Step 3 finds 


the sequence 3,(3,6),6,(6,8),8 as the Homing Path m, E 


E and s,(2,3),3 the Augmentation Path m = 
, , 


Step 6 sets to zero the branch flows а (3,2) from v, to v, and replaces 


3 t 
each branch flow in the path by zen. = 1, Step 7 truncates s'-3 path 


The path ew 1,3), 


Seas the Truncation Path m. 


by setting à (1,32) to zero, and Step 8 augments s-3 path by setting 
E 12,9) = Í. One s-t path is now processed and the algorithm will 
search foz another one until no more augmentation pathscan be found. By 
then, f will be maximal. 
SU 
. Flowchart of this al gorithm is in Appendix A-3 and its computer 


implementation MAXCON-1 is presented in Appendix B-1. 








3.3 ANALYSIS WITH RESPECT TO VERTEX DISCONNECTION 

The criterion of survivability attached to this problem is the number 
of stations that must be deleted from the network in order to disrupt com- 
munication. This is the smallest disconnecting set (0, of a network 
graph. Given any pair of vertices Ne and Vos assign unit capacities to 
all vertices except va and Me (which have infinite capacities), then the 


smallest s-t vertex cutset à is numerically equal to the maximum num- | 


‚t 


ber of vertex disjoint s-t paths. The value of б). t is the network ver- 
3 


tex-connectivity with respect to the given source and sink. 
3.3.1 s-t Vertex Connectivity 

To solve the problem of finding pe each"vervex іп G аи 
into two vertices connected by a branch of unit capacity forming a new 
graph G'. Algorithm 1 is then applied to G'. The result is a branch 
cutset composed only of the unit capacity branches. The regular network 
branches are | having infinite capacities. This approach is 
however very inefficient since the number of vertices are almost doubled, 

An algorithm known as Vertex-Pair Connectivity derived from Algo- 
rithm 1 was developed by Frisch ЕЛ which applies directly to graphs of 
directed, undirected, or mixed networks without splitting the vertices. 
The vertex-splitting process is absorbed in the modification of the label- 
ing routine of Algorithm 1. The algorithm is stated here in a different 
form, Several terms and rules must first be defined: 

1. A vertex is labeled with an ordered triple (i,+,M). The 
first and second entries are defined as in the original Labeling Routine. 
The third entry is a zero or one as will be explained in the algorithm. 

2. <A vertex is strongly labeled if the second entry is a minus, 
and weakly labeled if it is a plus. 

3. A circle around the second entry denotes the vertex is 


scanned, 
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4, If a vertex is weakly labeled and unscanned, and then 
receives a sirong label, the weak label is erased. However if the label 
is scanned, it should not be erased. 

Algorithm 4 

Labeling Routine 


1. Erase vertex labels. Label u- by (s,+,0). уа ls now 
labeled and unscanned. All other vertices are unlabeled and unscanned. 

2. Select any labeled "es vertex Vas lf none, ter- 
ninate. 

a. Suppose vi is weakly labeled: 

1) If £(1,3) = 0, (4,3) € f'uand de» o 
label all vertex Y; by (i,+,0). Y; is now weakly labeled and unscanned. 

2) If f(43) и ОН ВЕЕ а v; not 
strongly labeled, strongly label 77 by ae u is now strongly label- 
ed and unscanned, 

Scan Vi by encircling the + in the vi label. 

b, Suppose У. is strongly labeled 

1) if f(1,3)) = 0, (1,3) € Y ала Se. 
label all unlabeled vertices 2 by (i,+,M). у, is now weakly labeled and 
unscanned. | 

2) Tf (4,0) = 1, (CASA ЕЕ 2 not 
strongly labeled, strongly label Us by (1; 9M ER E is weakly labeled, 
erase the weak label, E is now strongly labeled and unscanned. 

Ih (DJ, if v, is both weakly and strongly labeled, 
set M = 1, otherwise M = 0, Scan vi by encircling the minus sign in the 
Vs label, 

3. If v, is labeled, proceed to Augmentation Routine. 


t 


Otherwise return to Step 2. 
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Augmentation Routine: 
ia bey z - "+ 
940 ТТ МЕ is weakly labeled by (q,+,M), set f(q,z) = 1. 
b. If v, is strongly labeled by (g,-,M), set Т 
3. If q = s return to Step i of Labeling Routine. Otherwise 
proceed to Step 4, | 
kk. Let z= = а. If M = 1 return to Step 2(b). Otherwise return 
Forsten 2(a). | 
To identify the s-t vertex cutset, let U be the set of all unlabeled 
vertices, and L be the set of all labeled vertices, Further let Г, be the 
set of labeled vertices connected to an unlabeled vertex by a branch di- 
rected from a labeled vertex to an unlabeled one, then the s-t vertex cut 
X is given by: 


а d L ee 


IM (L, - у.) U (nim с (ыш) Ta (О Ој 


Since the source vertex cannot be а component of the vertex cut, Ед. 


85.1.2) shows that if vo is connected to unlabeled vertices v,, then 


ү? 
the vertices v, are some component of the cut. 

The graph in Fig. 3.4 illustrates the application of the algorithm 
for a eel network. Assume that the path indicated by solid arrows 
is a result of initial labeling and augmentation. Assume further that 
the second application of labeling routine is just completed and Vy is 
labeled-as indicated by the vertex labels. When the augmentation routine 
is applied, it will backtrack the path in the usual manner, however, on 
reaching Vy it has to select the strong label since ve has M = 1. If the 


weak label were selected, a non-feasible flow would result. This explains 


the purpose of the third entry in the label. The only case where two 
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scanned labels can occur on a vertex is the situation in Ny where a flow 
of unit magnitude already exists through the vertex and this vertex re- 
ceives labels in the succeeding iteration of the algorithm. If the weak 
label in vj, Was unscanned when it received a strong label then the weak 
label should have been erased and МА would have been labeled (4,+,0). 

Program MAXCON-3A in Appendix B-6 is a computer implementation of 
Algorithm 4, A different way of circumventing the situation in vertex 
s of Fig. 3.4 was applied. For a vertex like Vy where a flow already 
exists, if in the succeeding labeling process Y receives a weak label, 
it will not be scanned, but if it receives a strong label it will be scan- 
ned. The latter would indicate that the existing outward flow from vj, can 
in effect be diverted to another path. In the former, it would be waste- 
ful to continue seeking another path passing through Vy when the vertex 
is already saturated. Іп no case then can a vertex acquire two scanned 
labels, The third entry in the label is therefore no longer necessary. 

The same form of label as in MAXCON-1A is used in this program. 
All steps in the algorithm that are associated with the third entry in 
the label are deleted. 

A flowchart of Algorithm 4 is in Appendix A-4, 
3.3.2 Multiterminal Vertex Connectivity 

The Flow Variation Algorithm of Frisch as presented in Sec. 3.2 
maximizes flows in a finitely branch-weighted network, and it was readily 
adapted to solving the multiterminal branch-connectivity. In this sec- 
tion the problem is different in the sense that the vertices control the 
information flow. 

In order to utilize the principle of the Flow Variation Algorithm 
here, the Routine A of the algorithm must be altered. Recall that Routine 


A is the Ford-Fulkerson Labeling Routine. If this routine is replaced by 
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the Labeling Routine of Algorithm + which is Frisch's Vertex-Pair 
Connectivity Algorithm, it would enable the Flow Variation Algorithm to 
handle the multiterminal vertex-connectivity problem. Consequently the 
method in which the s-g and s-t paths are augmented must also conform 
with the Augmentation Routine of Algorithm 4, 

With these changes, the algorithm flowchart is the same as in Ap- 
pendix A-3, however there is a significant change in the computer pro- 
gram MAXCON-3 as presented in Appendix B-5 Scien to MAXCON-1. 

The solution given earlier for obtaining the multiterminal vertex- 
connectivity оҒ а directed network applies equally well for undirected 
networks. ‘Sometimes one finds it impractical to compute all the entries 
in the entire connectivity matrix especially when investigating a very 
large network, It is often sufficient to just verify whether a network 
satisfies a certain minimum degree of connectivity. 

With this analysis criterion it is only necessary to verify that 
the network is at least R-connected. The connectivity being referred to 
here is the vertex-connectivity which we know is the smallest. 

The following theorems of Kleitman [14] greatly reduce the number 
of verifications necessary rather tnan doing A separate tests for 
the existence of R vertex-disjoint paths. 

Theorem: 

To verify that all n(n-1)/2 vertex pairs in G can be connected 
by at least R vertex-disjoint paths, it is only necessary to verify that 
some set of R vertices are each connected by at least R vertex-disjoint 
paths with every vertex in the graph, 

Two other theorems resulting from the above are stated here in an 


algorithm form and are in themselves the analysis procedure. 
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Algorithm 5 
1. Choose any vertex, say v, as source, 

Іеї і = 1, 

a. Select a vertex v. in Gi as sink, Remove from Gi all 
other vertices (v) ,k-71,...,L previously designated as sink that 
are adjacent to У. If 1 > ЕВ гео о Sup. 

1) For Ё connected to vi 
Excluding the direct en i-j, use Algorithm ^ to 
verify existence of q vertex-disjoint paths from Vi to RES Set s = an. 
If s«R go to Step !. Otherwise go to Step i(b), 

b, Are all n-i vertices processed? If not go to Step 1(a), 
otherwise go to Step 2. 

2. Decrease R өзін ЈЕ К = 0, terminate; otherwise go to 
Step 3. 


3. Remove Vi from G. e Choose vertex v. as source. Increase 


iti 
i by one, Return to Step 1(a). 
4, Record s-connectivity between Vi and ML Terminate.  Net- 
work is not R-connected, 
The first theorem may be restated in a more specific sense: 
If there are at least R vertex-disjoint paths between v4 and 
MI and v, and Ys in G, then to verify the existence of R such paths be- 
tween Vi and UE in G, we need only look for R - 1 such paths in G - {vi} . 
For the proof of these theorems and illustrative example, the read- 
er is referred to the original paper by Kleitman. 
А flowchart for computer implementation is given in Appendix А-5. 
Program MAXCON-4 in Appendix B-7 is encoded to analyze an undirected net- 


work for a uniform R-connectivity where the value of R is initially 


specified. 





IV. DESIGN OF LOW-COST SURVIVABLE NETWORK 


4,1 INTRODUCTION 

There are now a number of synthesis procedures for survivable networks 
appearing in literatures elsewhere. These are summarized by Frank and 
Frisch in [ 10 | and are categorized according to link or station damage. 
The synthesis approach are all aimed towards n NO of redundancy 
between stations and satisfying network redundancy requirement with the 
assumption that cost is either uniform for each link/station or not con- 
sidered at all. 

In the design of practical flow networks, the cost of construction 
increases with increase in redundancy therefore the design approach must 
take into account non-uniform cost for the construction of links or sta- 
tions, The goal is a minimum cost network that satisfies a prespecified 
redundancy requirement, So far there is no exact analytical method that 
could be applied to solve the problem. A procedure to search for a low- 
cost network with respect to vertex disconnection is presented by Steig- 
t2, etal [18] using an heuristic approach. The algorithms they devel- 
oped were based on the techniques for an approximate solution to the 
traveling salesman problem. The salesman's tour is a Hamiltonian circuit 


which is a special case of the flow network having a redundancy of two, 


4,2 DESIGN WITH BRANCH DISCONNECTION CRITERION 

In some cases branch destruction is more IHRE vertex destruc- 
tion and it is just as necessary to obtain a design method using the 
former criterion as the latter. A procedure is presented in this chapter 
for design of a low cost network (with high probability of being optimal) 


with at least г, | branch disjoint paths. 
) 
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The same procedure used by Steiglitz, etal is employed consisting of 
two parts. Тһе first part called the Starting Routine generates a fea- 
sible initial network (one which satisfies redundancy requirement) and 
the second part is the Optimizing Routine which makes an iterative im- 
provement on the network cost by branch exchanges while maintaining fea- 
sibility. Portions of the procedure that differ from that of Steiglitz 
are the feasibility testing of the network and the optimizing routine. 
The method of optimization that will be presented is an adaptation of the 
recently published heuristic algorithm for the traveling salesman problem 
by Lin and Kernighan [16]. 

The heuristic design algorithm is stated as follows: 

Algorithm 6 

1. Generate a feasible initial network G that satisfies the 
specified redundancy R using the Starting Routine, 
2. Attempt to find an improved feasible solution G' by some 
local transformation using the Optimizing Routine. 
3, If a similar network with a lower cost is found, then re- 
place G by G' and repeat from Step 2. 
lh, If no local improvement can be found, G is a locally optimal 
solution. Repeat from Step i until computation time runs out or the so- 
ita on is satisfactory. 
4.2.1 The Starting Routine 
This routine was developed by Steiglitz, etal in jubes) 
1. Plot vertices representing stations in a Euclidian 2-space. 
2. Number the vertices uniformly at random. 
3. Assume cost to be proportional to distance and obtain the 
cost matrix from the Euclidian distances between vertices truncated to an 


integer. 





!. Assign redundancy to each vertex, 

5. Add branches to the network one at a time. Reduce by unity 
the assigned vertex requirement every time a link is connected to that 
particular vertex, A branch is added in this manner: 

a. Connect the vertex with highest updated requirement to 
another vertex with the next highest updated requirement; for all vertices 
having equal requirement, choose one that results in lower cost. 

b. In cases of ties in ET and cost, these are 
resolved by choosing the vertex highest on the list. 

с. No parallel branches are allowed. 

6. When the requirement in all the vertices are reduced to 
zero, test feasibility of the initial solution using Algorithm 2. 

7. If solution is feasible go to Optimizing Routine, otherwise 
repeat from Step 2. 

This routine is made nondeterministic by numbering the vertices 
uniformly at random at each start. Randomized starting solution has been 
observed experimentally to produce a variety of starting networks as con- 
trasted to constructive solutions which are deterministic where more than 
one solution may not be obtained, 

-The analysis tools presented in Chapter III come in handy for test- 
ing feasibility. For the initial network, the one that fits is Algorithm 
2 which determines the overall branch connectivity of the network. А 
complete checkout need not be carried any further once infeasibility is 
detected. 

4.2.2 Optimizing Routine 

The optimizing procedure applied by Steiglitz is the X-change where 

two branches are interchanged at a time, and determined whether the trans- 


formation is favorable and preserves feasibility. Lin and Kernighan have 
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lately experimented on a more systematic and effective way of branch 
transformation in a network with 2-redundancy. Their results are applied 
here for a network with redundancy greater than two. 

The basic ideas underlying the transformation are as follows; 
Consider a nonoptimal but feasible network G with cost c(G) and an opti- 


mal or near optimal network G' with the same number of vertices and with 


cost c(G') < c(G). G is nonoptimal because it has К links Xen, that 
are "out of place". TEE G optimal they should be replaced by k links 
Ууз... зуг of G'. We try to find k by identifying sequentially, element 

by element, two disjoint sets of links X - eon] and Y = E 


using a gain criterion. Let the costs of X. and Yi be |х, | апа | xil 


respectively. Define g, = х, | - |у. | as the gain from exchanging x. 
with У; • For some i, E may be negative however we want the final gain 
to be positive, i.e., > б, > 0. We M looking for sequences of 6,5 whose 
partial sum is always ке ы When 2 & X О, the search is stopped. 

i 


The k links in X are deleted or "broken" and replaced by k links 
in Y resulting in G' which has lower cost. The process is iterated until 
no further reduction in cost is attained. 

We now state the routine more formally assuming that a starting 
solution has been generated: 

1. -Set P* 2 0, (P* is the best improvement so far). Choose 
any vertex as Uu and let Xy be one of the edges of G incident to tye Let 
EE 1, 

2. Designate the other endpoint of x, as tos From t choose 
y, different from existing links in G such that g,>0 and designate the 
other endpoint of y, as tae If no such МИ exists, go to Step 5d). 

3. Leti =i + 1. Choose X. (which currently joins tos 1 to 
+.) апа y, as follows: 
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а. х is chosen so that it is the longest link іп G 


incident to vertex t.. 
21-1 


b. Y. is some link that can be drawn from vertex toi sub- 


ject to (c), (d), and (e) below. If no y, exists, go to Step 4, 
| c. To guarantee that x's and y's are disjoint, x, cannot 
be a link previously joined (i.e,, a y j « i), and similarly y, cannot 
be a link previously broken, 
д. Р, = z 8, > 0. (Gain criterion) 
е. Уу, chosen must permit breaking of an Х 44 


f. Before МА is drawn, check if closing up by joining 


tos to t, will give a gain value better than the best obtained previously. 


E А А ¥ = = х 
Let y,* be a link connecting t,, to t, and let g, х, | |У; | 4 


х х х = 
EE + а “> р, деррж— Ре 


best improvement in G. P* » O and must be monotonically nondecreasing, 


+ g.* and let k =i. (P* is now the 
ل‎ 


The index k defines the sets to be exchanged to achieve P*). 


4, Terminate the search for X. and У; in Steps 1 to 3 when 


either 
a. no further links x, and y, satisfy 3(c) to (e). 
DD P 
Pc 
Delete the x, 's and replace them by the У;'5, ЕЕ: 
tO obtain G'. 


5, Test the feasibility of G' (see Sec., 22200, ШО 7 
take G' as the starting network and repeat from Step 1. Otherwise gen- 
erate another initial solution using the Starting Routine, 

6. If P* = 0, a limited backtracking is invoked: 

a. Repeat Steps 3 and 4, choosing Y'S in increasing cost 


as long as they satisfy the gain criterion 54 t g, > 0. 
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b. If all choices of Yo in Step 3(b) are exhausted without 
profit, return to Step 3(a) and try the alternate choice for X2 
с. If this also fails to give improvement, " further back- 
up is performed to Step 2, where the y,'s are examined in increasing cost. 
d. If the y¿'S are also exhausted without profit, try the 
alternate Xy in Step 1. 
e, If this fails, select another vertex ty and repeat 
from Step 1. 
7. The routine terminates when all vertices in G have ән 
designated ty without profit. Return to the Starting Routine, 
A justification of Step 6 is in order: Experimental results 
by Lin and Kernighan showed that for the 2-redundancy network, backtrack- 
ing on the third level and higher, i.e., i > 3, indicated considerable 


time penalty, hence the only alternates for х1, Уң апа х У, at the first 


2 
and second levels are examined, Other refinements to the procedure dis- 
cussed in the original paper are not applied here since actual computer 
runs may indicate different time-cost trade-offs for networks with redun- 
dancy greater than 2. 
4.2.3 Feasibility Test 

“The optimizing routine for 2-redundancy network has a built-in 
feasibility test but is not applicable to higher redundancy networks, 
The fact that G is feasible does not imply that G' is also feasible even 
though the degree of the vertices satisfies the requirement. It is there- 
fore necessary to test the redundancy after every transformation. Rather 
than applying Algorithm 2 to examine the new network G', the number of 
flow calculations can be reduced significantly using the theorems stated 


by Steiglitz, etal for an X-change or 2-change. Extending this to a k- 





change transformation and with uniform redundancy, we need only to check 
the redundancy between the pairs of vertices whose connecting links were 
deleted, Algorithm 1 is most suitable for investigating feasibility after 


every any transformation. 


Example: 

Given a set of stations whose geographical arrangement is trans- 
ferred on an Euclidian plane as in Fig. 5.1. We attempt to construct a 
minimum cost flow network with a uniform redundancy of 3. We model the 
network as a graph of vertices and branches and illustrate the applica- 
tion of the design procedure, 

otart by numbering the vertices uniformly at random as in Fig. 5.1. 
Measure the distances between vertices (assuming proportional to cost) and 
obtain the cost matrix. Add branches one at a time and keep track of the 
number of incident branches on each vertex by using an updated require- 
ment table as shown. Row 1 of this table is filled up by the assigned 
vertex requirements. This row has all the same entries so start with 
vertex 1. The cost matrix shows branch (1,6) has the least cost (here 
the parentheses are used in place of brackets to represent undirected 
branches).  Encircle T in the matrix and construct the branch 
M 6) in Fig. 5.2. Form row 2 of the requirement table by reducing en- 
pres for mic 1 and 6 by unity. The highest updated requirement in 
this row is 3 and the highest on the list is vertex 2. From the cost ma- 
trix we see that branch (2,7) has the least cost.  Encircle this entry 
and construct the branch (2,7). The procedure is continued and the se- 
quence of branches constructed next are (3,4), (5,8),(1,7),(2,8),(3,6), 


(4,5) 2) 2:7) , (4,8) , and (5,6) e 
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When the number of vertices is odd, the last row of the requirement 
table has one vertex left with an updated requirement of 1 while the rest 
are zero. Connect this vertex to one which produces the least branch cost, 
at the same time satisfying Step 5(c) of the starting routine. 

The initial solution generated by the starting routine is shown by 
solid lines in Fig. 5.2. It is found feasible using Algorithm 2. The 
cost of the network is obtained by summing up the values of the encircled 
entries in the cost matrix. For this A the initial созба СІ 

We apply now the Optimizing Routine to reduce the cost. Set P* = 
О ала 1 = 1. In Fig. 5.2, choose vertex 1 as ti’ x is selected such 


that it is the longest among the three incident branches to tye Vertex 


2 which is the other endpoint of Xy is designated t Brom bg RU Y4 


2° 2 

using Step 3 of the routine such that 84 > 0. Only the imaginary branch 
(2,4) is available since the other possibilities (2,5),(2,6), and (2,3) 
would result to в, < 0. Choosing (2,4) as y, the gain g, = ES - 


|у; | = 3, The overall profit Р, = 3 > P*, Designate vertex 4 as Las 
Now let i = 2, Choose X, as indicated.  Assign ty to vertex 8 


which is the endpoint of X». Using Step 3 again, Yo is selected to be 


the imaginary branch (8,7). The gain 6 = 9 апа Р, == Р, 


fore going further, check if joining ty to ti would improve P*, With 


T 65 = 12, Be- 


(8,1) as Y5*, Ea — 9 and P* - Р, t 8,” = 12. We have to terminate the 


‚search for x,'s and y,'s since P* =P, (Step 4). Had P* been greater 


2 
for X4 and Уа” 


than P, we could have continued by designating vertex 7 as Ue and look 


Construct the Y's and delete the X's, The network with improved 
cost of 249 (not yet optimal) is shown in Fig. 5.3. The cost is obtained 
by subtracting the costs of X, and X, from the initial value and then 
adding the costs of the branches У; апа Y5*. 
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Feasibility of the new network is tested using Algorithm 1 since 
we need only to check the branch connectivity between the vertex-pairs 
1-2 and 4-8. This network is feasible so we can iterate the procedure 
using Fig. 5.3 as the new starting solution, 

Satisfying the vertex requirement does not necessarily mean fea- 
sibility is preserved. Fig. 5.4 obtained by the optimizing procedure is 


not feasible as it has a minimum connectivity of 2. 
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V. CONCLUSION 


Analysis procedures provide convenient tools for evaluating the 
"strength" or "weakness" of flow networks with respect to some physically 
meaningful criteria. The computer analysis programs were prepared as im- 
plementation of the algorithms developed and are capable of investigating 
large networks based on either branch or vertex survivability measure. 

It was shown that the analysis tools are not only used for measuring 
performance of existing systems but are also applied effectively in their 
design particularly when either explicit solutions are impractical or so- 
lutions do not appear to be obtainable. A heuristic algorithm for the 
design of minimum-cost survivable networks is presented. Analysis tech- 
niques are found indispensable in the procedure. 

Further study is indicated in design perhaps by encoding the algorithm 
for computer running and to ере the following: 

1. What level of backtracking would be considered optimal? 

2. How often does infeasibility occur in the starting routine? The 
Peovability of occurrence "e be different from that obtained by Steiglitz 
since the criteria are not the same. 

3. Are refinements described іп [16 | applicable to networks with 


redundancy greater than 2? 
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APPENDIX A 


FLOWCHARTS 


s-t BRANCH CONNECTIVITY 
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Àny more 
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А-2. MULTITERMINAL BRANCH CONNECTIVITY 


( Undirected Network ) 


Read Connection 
Matrix 


lnitialize vari- 
ables matrices, 
source, and sink | 






Find initial 
s-t cut, 
Use Algorithm i 
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List down labeled 
vertices under X and 
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XBAR has single ИЕСИ X has single ЊУ 
vertices are in 
X and XBAR 
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Condense XBAR, 
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Find s-t cut 

. Use Algorithm 1 
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List down single 
vertex in X as 
processed; 










List down single 
vertex in XBAR 
as processed, 








A-2, (continued) 











Are all 
XBAR vertex-grps 
processed? 
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А-3. MULTITERMINAL BRANCH CONNECTIVITY 


( Directed Network ) 
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A-4, s-t VERTEX CONNECTIVITY 
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А-5. UNIFORM R- CONNECTIVITY 


Read, N, R, 
connection matrix 


Choose source 
vertex - 
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